Relativistic Weierstrass random walks 
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The Weierstrass random walk is a paradigmatic Markov chain giving rise to a Levy-type superdif- 
fusive behavior. It is well known that Special Relativity prevents the arbitrarily high velocities 
necessary to establish a superdiffusive behavior in any process occurring in Minkowski spacetime, 
implying, in particular, that any relativistic Markov chain describing spacetime phenomena must be 
essentially Gaussian. Here, we introduce a simple relativistic extension of the Weierstrass random 
walk and show that there must exist a transition time tc delimiting two qualitative distinct dynam- 
ical regimes: the (non-relativistic) superdiffusive Levy flights, for t < tc, and the usual (relativistic) 
Gaussian diffusion, for t > tc. Implications of this crossover between different diffusion regimes are 
discussed for some explicit examples. The study of such an explicit and simple Markov chain can 
shed some light on several results obtained in much more involved contexts. 

PACS numbers: 05.40.Fb, 47.75. +f, 95.30.Lz 



I. INTRODUCTION 

Relativistic Brownian motion is an interesting and ac- 
tive area of research nowadays [1|. The dynamical be- 
havior of any relativistic process is, of course, strongly 
constrained by the main properties of the relativistic 
kinematics, namely the speed of light c as the maximum 
possible physical velocity, the invariance under Lorentz 
transformations, and the causal structure associate with 
the light-cone. Some well established non-relativistic dy- 
namical behaviors are simply incompatible with the prin- 
ciple of Special Relativity. This is specifically the case 
of the so-called Levy flights [2|, where very rare events 
with arbitrarily high velocities give rise to a superdiffu- 
sive regime characterized by a power law 



(^'(i))«i^ 



(1) 



where /it > 1 is the corresponding anomalous diffusion ex- 
ponent. However, Special Relativity does not allow such 
arbitrarily high velocities, preventing the appearance of 
such superdiffusive regime in Markov chains involving 
spacetime events. These points are studied, for instance, 
in [3| by using relativistic versions of the Fokker-Planck 
equation and of the fluctuation-dissipation theorem, lead- 
ing to a fractional-derivative extension of the diffusion 
equation. An earlier analysis of relativistic random walks 
can be found in [J]. In Q, a generalized Wiener pro- 
cess avoiding superluminal propagation is introduced, 
giving rise to a non-Markovian relativistic diffusion pro- 
cess. The influence of the spacetime causal structure on 
dynamical processes, namely the implication of the pres- 
ence of an event horizon, was investigated recently in [6| . 
We consider here a simple relativistic extension of the 
Weierstrass random walk in order to shed some light 
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on several relativistic aspects of Levy flights. We re- 
mind that the usual one-dimensional Weierstrass random 
walk[7j corresponds to a Markov chain governed by the 
following probability density function 
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with the "jump" function 
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where a > 1 and 6 > 1 are dimensionless constants and 
V > gives the scale of the jumps. A particle moving ac- 
cording to (m and ([3]) can perform jumps, in both direc- 
tions, with magnitude v,bv,b'^v,b^v, . . ., and with prob- 
ability, respectively, given by (a — l)/a, (a — l)/a^, (a — 

l)/a^, (a — l)/a'^, We consider here that each step i 

lasts for a given and fixed time interval and, hence, the 
asymptotic dynamics for large times and for large num- 
ber of steps are identical. One can think the constant v 
as the velocity acquired by the particle, just prior to the 
jump, by some unspecific microscopic mechanism. No- 
tice that one can indeed define a continuous-time version 
of the Weierstrass random walkfq, but for our purposes 
here, the simple Markov chain governed ^ is enough. 

We are mainly interested in the anomalous diffusion 
process associated to the Weierstrass random walk. A 
brief review of the main results concerning this topic is 
necessary here. We wish to characterize the large L be- 
havior of (Ar^(L)), where 



X{L) = Xi + X2+ X^ -^ ■ ■ ■ -{- XL, 



(4) 



with X( being the size of the l*"^ jump, with probability 
density function given by ([2]) and (|3]). Since xi are inde- 
pendent random variables for different i, we will have 
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Notice that 



X 4'{x) dx 
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from where we have that, for b^ < a, {x'^) is finite, equa- 
tion ([5]) reduces to the usual diffusion process with /x = f , 
and, thanks to the central limit theorem, the total proba- 
bility density function describing the passage from £ ~ 1 
to £ — L will be very close to a Gaussian. For b^ > a, 
on the other hand, {x"^) diverges, the central limit theo- 
rem cannot by applied anymore and the associated dis- 
tributions will not be Gaussian, leading to an anomalous 
diffusion exponent /i > 1. We can use the heuristic ap- 
proach of [3 in order to evaluate fj, for this case as well, 
which correspond to a < 2, where 



a 



logQ 
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The key idea is that, according to ([2]) and ([3]), if one 
considers a large number L of steps, the most frequent 
steps will have magnitude v. Steps with magnitude bv 
will be 1/a less frequent than the previous one. In gen- 
eral, steps with magnitude fe"+^ti will be 1/a less frequent 
than those ones with magnitude 6"u. Also, if we consider 
large, but finite, number of steps, the summation ([2]) will 
be effectively truncated at a given value n = n^ax- Let us 
consider the following succession of steps, which exhibits 
the required hierarchy of jumps. 
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1, 
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corresponding to a"'"^"' steps with magnitude v, a""""'"^ 
with with magnitude bv, and so on. For large nmax, we 
have 



3=0 
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Yet for large but finite nmax, we can estimate (a;^) as 



(x2) _^ a_i"— /62XJ 
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For a > 2, as we already know, we have the Gaussian 
result since 



(x^) 
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which coincides with the exacted result evaluated from 
^. For a < 2, on the other hand, equation (|10p implies 
that 
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where (0) was used, leading to a superdiffusive {fi — 
2/ a > 1) behavior characterized by 



(X^L)) /a-1 
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The relativistic kinematics, however, will change dramat- 
ically this scenario. 



II. THE RELATIVISTIC WALK 

The microscopical origin of the jumps in a non- 
relativistic random walk is not relevant from the dynam- 
ical point of view. Provided that the Markov property 
holds, i.e., the position X{£) of the system at a given 
step £ depends only on the position at the previous step 
X{i~l), and that successive jumps are independent ran- 
dom variables, the standard approach to lead with ran- 
dom walks [1| can be applied. Let us, nevertheless, sup- 
pose that the jumps are due to microscopical collisions as 
in the historical example of Brownian motion. The hier- 
archy of jumps V, bv, b'^v, b^v, ... is, in the nonrelativistic 
case, associated with a similar hierarchy of momentum 
transfers in the collisions p, bp, b^p, b'^p, . . ., withp — mQV, 
where mo is the particle rest mass. However, according 
to Special Relativity, the velocity and the momentum of 
a particle should obey 



P = 



niov 
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where /3 — v/c, with profound implications for the jump 
hierarchy. Assuming that the same hierarchy of momen- 
tum transfers is present in the relativistic case, we will 
have the following jump function 



^n — 'J n 



6" 



y/TTWi^' 



(15) 



instead of ^. The hierarchy of relativistic jumps will 
be {Jnv}, occurring in the dynamics, respectively, with 
probability {(a - l)/a"+i}, with n = 0,1,2,.... No- 
tice that, in frank contrast with the nonrelativistic case, 
for large n, one has JnV ~ c, meaning that there will 
be no arbitrarily large jumps in the relativistic case, in 
agreement with the fact that no acquired velocity by any 
microscopical mechanism can exceed c. The first conclu- 
sion we can draw from this relativistic extension of the 
Weierstrass random walk is that the whole process must 
be Gaussian, since, in this case, we have 
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for any /3 > 0, implying the usual diffusion with /i = 1, 
irrespective of the value of a. 

Typical non-relativistic situations are characterized by 
a small /3. For such cases, from (fT6| . one realizes that 
there should exist a critical value ric such that, for 
n ^ Tic, the summand of ([T5|) is essentially the same 
one of the nonrelativistic case. If it is possible to choose 
a large n^ax "^ n^, one could apply the same heuristic 
approach of last section, implying in a superdiffusive be- 
havior with /I = 2/a > 1. (We assume hereafter that 




FIG. 1: Diffusion in the relativistic Weierstrass random walk. 
The sohd line (left scale) corresponds to (|19|l . which, for large 
rimax, is well approximated by (|20[l . The dotted line (right 
scale) is the ratio between the non-relativistic ((Sjl and rel- 
ativistic pSjl jump functions. It is clear the appearance of 
a crossover between the different diffusion regimes near Lc 
given by (|18|) . For L < Lc, we have Levy- type superdiffusive 
behavior, while for L > Lc the dynamics settle into a Gaus- 
sian diffusion. This plot corresponds to the case where a = 2, 
b = 3, and^= 10""'. 



< a < 2.) This occurs because L is large enough to 
justify the averages of last section, but it is still small 
enough to guarantee that the large n events that would 
imply the relativistic regime are extremely rare and will 
not contribute effectively to the averages. In other words, 
the system needs some time to realize that it is indeed 
relativistic! If we allow for rimax 3> Uc, we will have in 



((161) the convergent relativistic summation, implying the 
usual diffusion with exponent fi = 1. It is clear that 
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where ([9|) was used. If /3 is small, as one expects in the 
typical nonrelativistic problems, Lc will be large. For 
L <^ Lc, the system behaves as in the nonrelativistic 
regime and exhibits the properties of a Levy flight with 
/i = 2/a. On the other hand, for L 3> Lc, the system 
change its behavior to the relativistic regime, character- 
ized by ordinary diffusion with /i = 1. Such crossover 
between different diffusion regimes, depicted in Fig. [1] 
is compatible with the results of ilC)]. In fact, our model 
can be viewed as a simple microscopical realization for 
the generalized Fokker-Planck equation considered there. 
We can apply the same heuristic approach of last sec- 
tion to the relativistic case. In particular, we have for a 
finite and large rimax 
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leading to the following generalization of ([T3| for /3 7^ 
(see the Appendix for details) 



(X^L)) 
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where F(a, b; c; z) stands for the standard hypergeomet- 
ric function[ll[. Using that F{a,b;c;0) — 1, we have 
from (|20p an anomalous diffusion process with ^ = 2/a. 
for large L obeying L <C Lc- On the other hand, for 
L » Lc, we have (see the Appendix for details) 



(21) 




leading to the usual Gaussian diffusion 



(XHL)) 
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for L:^ Lc 



III. DISCUSSION 



The results of the preceding sections can be summa- 
rized as follows. Suppose we have a Weierstrass ran- 
dom walk model with typical velocity v, implying in Levy 
flights characterized by an anomalous diffusion exponent 
/i > 1. Then, relativistic effects imply that, after a cer- 
tain critical number of steps Lc ~ (c^/u^)", the system 
loses it anomalous diffusion properties and the dynam- 
ics necessarily settle into a Gaussian diffusion. In order 
to estimate the order of magnitude of these relativistic 
effects, let us associate the kinetic energy of the walk- 
ing particle, and thus v, with the typical thermal energy 
ksT . For non relativistic situations, where fc^T is small 
if compared with ra^c^ , we will have niQV^ /2 w fc^T, 



leading to 



Lr 






(23) 



where variable w = 6^" was introduced and a is, in gen- 
eral, an irrational. For \l3^w\ < 1, we can introduce the 
following series expansion 



As our first explicit example, let us consider a system 
composed by helium atoms, for which ttiqc^ /ks ~ 4.36 x 
lO^^K. For such a system at room temperature (T w 
300K), a ballistic (/i = 2) Levy flight originated in a 
Weierstrass random walk, will become Gaussian due to 
relativistic effects after Lc ~ 2.7 x 10^ steps. Helium 
atoms at the surface of the sun (T w 5 x lO'^K) can 
experiment ballistic Levy flights in a Weierstrass random 
walk for no more than Lc ~ 6.6 x 10'* steps. In the interior 
of the sun (T « 5 x lO^K), only L^ ~ 2,000 steps wiU 
be enough for the the dynamics settle into a essentially 
Gaussian regime. 

Heavier particles or bodies will, naturally, lead to 
larger values for Lc- Let us take, for instance, the case 
of an Escherichia coli bacterium, for which toq — 665 
femtograms T^, leading to rriQC^/kB « 4.33 x lO^^'K. At 
room temperature, Lc « 4.1 x 10^^ for an Escherichia 
coli undergoing a Levy flight originated in a Weierstrass 
random walk with diffusion exponent close to those ones 
observed experimentallv|13| in systems of breakable mi- 
celles (/x K, 1.4). One realizes that such microscopic bod- 
ies can indeed experience much longer anomalous diffu- 
sion process than atomic scale particles. 
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and the integral (jAip will be written as 
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where F{a, b; c; z) is the standard hypergeometric func- 
tion [ll| . Notice that the hypergeometric functions have 
a single valued analytical extension over the entire com- 
plex plane, with the only exception of the positive real 
axis for z > 1 [11| , justifying the use of (|A3[) for the eval- 
uation of the integral (jAll) . which limits, in fact, do not 
belong to the region where the expansion (jA2p converges. 
The integral (lAip may also be evaluated by exploring the 
hypergeometric function identities [ll| 
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Appendix A 

We can approximate (1191) by an integral by using 
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and F{1, b; b; z) = (1 — z) ^, valid for any b and c. Equa- 
tion (PH)) follows straightforwardly from (jA3[) . 



The evaluation of ([20)1 for L ^ Lc requires an asymp- 
totic analysis for the hypergeometric function. By using, 
for instance, the identity 15.3.8 of [ll[, we have 
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implying that, for large z and < a < 2, 
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leading finally to (HI]). 
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